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Abstract 

We compute the large scale (macroscopic) correlations in ensembles of normal 
random matrices with a general non-Gaussian measure and in ensembles of general 
non-Hermition matrices with a class of non-Gaussian measures. In both cases 
the eigenvalues are complex and in the large limit they occupy a domain in the 
complex plane. For the case when the support of eigenvalues is a connected compact 
domain, we compute two-, three- and four-point connected correlation functions in 
the first non- vanishing order in in a manner that the algorithm of computing 
higher correlations becomes clear. The correlation functions are expressed through 
the solution of the Dirichlet boundary problem in the domain complementary to 
the support of eigenvalues. 
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1 Introduction 



A matrix M is called normal if it commutes with its Hermitian conjugated: [M, M^] = 0. 
The both matrices can be simultaneously diagonalized, the eigenvalues being complex 
numbers. The partition function of normal matrices has the general form 

Zn=( d/i(M)e^*'^^(^^'^'). (1) 

J normal 

Here % is a. parameter, and the measure of integration over normal N x N matrices is 
induced from the flat metric on the space of all complex matrices. Like in other random 
matrix models, the large limit of interest implies oo, h 0, while A^^ stays 

finite. 

In a particular case, when Laplacian of the potential W{z, z) is a constant in a big 
domain of a complex plane, i.e., 

W = -MM^ + V{M) + V{M^), (2) 

where V{z) and V{z) are holomorphic functions, the normal matrix ensemble is equivalent 
to the ensemble of general complex matrices. It generalizes the Ginibre-Girko Gaussian 
ensemble . In this, perhaps the most interesting case for applications, the model bears 
some formal similarities with the model of two Hermitian random matrices j2] and the 
matrix quantum mechanics in a singlet sector [31 Unlike models of few Hermitian 
matrices, the normal matrix model is integrable for a general class of potentials, not only 
of the form (j21). 

Applications and studies of matrix ensembles with complex eigenvalues are numerous. 
A large list of references can be found in recent papers 0. New applications to diffusion 
limited growth models (Laplacian growth) [H], complex analysis |Zl-[ini and Quantum 
Hall effect jTT] were found recently. 

Despite of a comprehensive literature on the model of one and two Hermitian matrices, 
interest to the normal matrix ensemble, first introduced in and further studied in 
[13 , as well as to a model of general complex matrices with potential just starts to 
grow. In this paper we revisit the ensemble of normal random matrices to calculate the 
large A^ limit of correlation functions under condition that separation between arguments 
is much more than an average distance between eigenvalues (macroscopic, or smoothed 
correlations). Short scale (microscopic) correlations in the ensemble with the potential 
(j2)) are well studied (see e.g. for a review). 

At large A^ the eigenvalues of the random matrices are distributed within a domain 
(with sharp edges) of a complex plane with a density proportional to Aiy, where A is 
the 2D Laplace operator. We assume that the support of eigenvalues is a single connected 
bounded domain and that the boundary is a Jordan curve. Analytical properties of 
this curve determines the correlation functions. 

We will show that correlation functions are expressed through the objects of the 
Dirichlet boundary problem of the domain complimentary to the support of eigenvalues. 
Namely, the two-point correlations are expressed through the Dirichlet Green function, 
while higher multi-point correlations are expressed through the Neumann jump on the 
boundary, through the Bergman kernel and through the curvature of the boundary. The 
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objects similar to the correlation functions previously appeared in studies of thermal 
fluctuations in classical confined Coulomb plasma fT^ (see also pE| for a review) and 
in recent studies of integrable structure of the Dirichlet boundary problem in If 
the potential W is such that the support of eigenvalues collapses to a cut (or cuts), the 
normal matrix model reproduces known large A^-limit features of the Hermitian matrix 
model. 

We compute two-, three- and four-point connected functions of density of eigenvalues 
in the leading order in N, in a manner that the algorithm of computing higher correlation 
functions becomes clear. Connected density correlators are localized at the edge of the 
support of eigenvalues and show some universal features. They depend on the potential 
only through a shape of the support of the eigenvalues and boundary values of a finite 
number of derivatives of the potential. In the literature on Hermitian random matrices 
the universal character of correlations has been emphasized in Refs. [THl I17j. The two- 
point functions are distinguished since they depends on the support of the eigenvalues 
only. The two-point function has been previously computed in Refs. [Hj in the course 
of studies of confined Coulomb plasma. 

Other aspects, including integrable structure of the model (discussed already in [T31 
ITHj). semi classical properties of biorthogonal polynomials, critical points (a boundary 
with cusps), disconnected supports of eigenvalues, correlations at short distances, etc., 
will be addressed elsewhere. 



2 Preliminaries 

The measure of normal matrices. The measure in (P) is induced from the flat metric 
||5M|p = tr {SMSM"^) in the space of all complex matrices. Formally, one can write 
dfi{M) = 6{[M, M'^])dMdM'^ , where the delta-function selects the subspace of normal 
matrices. To introduce coordinates in this subspace, one uses of the decomposition 
M = UZU"^ , where f/ is a unitary matrix and Z = diag [zi, . . . , zn) is a diagonal matrix 
of eigenvalues . The measure is then given by 

*w-i^wk)'^"<^>iin<'^^.. (3) 

here d'^z = dxdy for z = x + iy, dfiQ is the Haar measure on the unitary group U{N), 
and An{z) = det(2;*~^)i<jj<Ar = Y\fyj{zi — Zj) is the Vandermonde determinant. 

If A{M) is any invariant function (i.e., a symmetric function of eigenvalues) of a 
matrix M (and M^), then the mean value 



Jdfi{M)A{M)e 



(A) = ' ' ^ / (4) 



is expressed through the integral over eigenvalues 

N 

■ I A{z)\A^{z)\' 



{A) = ^ / A{z)\A^{z)\' n (e^^(^-^"^-)ci%-) , (5) 
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where the partition function is 

1 r ^ 

ZN = j^J\^N{z)?\{{e^''^'^-^^d^z,). (6) 

Potential. For notational simphcity, we shall write simply W{z) instead of W{z,z). 
We assume that is a real-valued function with (at least local) minimum at the origin 
and set W{0) = for convenience. We also assume that the integral (jH)) converges and 
that is a regular function in both variables at the origin. We shall also set 

a{z) = --dAW{z) = -^AW, (7) 

and assume that cr(z) > 0. 

A special interesting case JT] occurs if the potential is harmonic in some big domain 
around the origin. Then in this domain a{z) = constant, say, set to be vr"^, and 

W = -Izl"^ + V{z) + V{z), (8) 

where V{z) is a holomorphic function. 

The measure of complex matrices. If the potential is chosen in the form (jSJ, then 
the same measure, up to a numerical factor, appears in complex random matrices. In 
this case the relevant decomposition reads: M = U{Z + R)U\ where U , Z are again 
unitary and diagonal matrix respectively, and R is an upper triangular matrix. The 
measure (jH} acquires a multiplicative factor Hij dRije~^^J (we use the fact that tr MM''' = 
XxZZ'^ + \xRR)\ As a result the representation (0) holds |19j . 

Coherent states of particles in magnetic field. Most of the known matrix ensem- 
bles represent coherent states of fermions. Fermonic representations are well known 
in the case of Hermitian matrices (see, e.g., [20]); where fermions live on a line and are 
confined by a potential. 

Complex matrix ensembles also enjoy a fermionic representation j^UE]- In this case 
fermions are situated on a plane and occupy the lowest energy level of a strong magnetic 
field B{z), the magnetic field is not necessarily uniform. It is related to the potential by 
B{z) = 2t[(j{z). The coherent state of the fully occupied lowest level is 

^{z^,...,zj,) = ^AN{z)e^U^-'^^'^\ Z^= \\^{z^,...,z^)\^^d\,. (9) 

So Zjv, in this picture, is the normalization factor of the wave function. 

Example: Gaussian model. In the case 

W{z) = -TX(y\zf ^21le{txz + t2Z^) , a > , 
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(the Ginibre-Girko ensemble) the partition function can be found exphcitly even for finite 
(see e.g. j2I]): 

ZN{a,ti,t2) = Zy{na)'2^ '{n a - At2t2) 2 exp — — — ^ , (10) 

\ h TX^a^ - At2t2 ) 

where = h^^^^^^^Tc^ Y[f=i j^- is the partition function for W = —\z\^. Correlation 
functions in this case are expressed through Hermite polynomials and are known explicitly 
for any A^. 



Correlation functions of traces. In this paper, we are interested in correlation func- 
tions of products of n traces, (nr=i fii^))- Clearly, they are expressed through ra-point 
correlation functions of the density of eigenvalues, 

N 

1=1 

We obviously have: 

(n \ „ n 

ntr/,(M))= / (p(^i)---P(^n))/l(^l).../n(^n)n A- (H) 
i=l I j=l 

So, the density correlation functions carry the necessary information. 

It is customary to deal with connected part of a correlation function. In case of the 
2-point function it is 

(p(2l)p(22))conn = (p(^l)p(^2)) ' {p{zi)) {p{z2)) . 

As N 00, the n-point correlation function of densities is 0(1) while the connected 
part of the n-point function is 0(iV^~^"). 



The method of functional derivatives. Correlation functions are the linear re- 
sponses to a small variation of the potential. Variation of the partition function (jUj) 
over a general potential W{z) inserts ^ S^'^\z — Zi) into the integral. Then 



_2 6\og Zn I , . , ^^ ^2^ (P(^i)) 



("W'^^' W- W->)Pfe))„„„ = s^^^ . (12) 



and, more generally, 

(p(Zl) . . . p{Zn))^^^^ = h 



2n log Zn 

6W{zi)...6W{z„ 



We shall use this method in the following version. Set 6W{z) = eg{z), where g is an 
arbitrary smooth function on the plane and e 0. Then, in the first order in e, 

h'5{tTfiM))=e{trfiM)tTgiM))^^^^. (13) 
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It is often convenient to consider correlations of the potential 

ip{z) = J log \z - z'\^p{z')dh' = htT {\og{z ~M){z- M^)) , (14) 
and of the current field 

J{z)^d^{z)=htr[^-^y (15) 

rather than correlations of density. The potential is the Bose field or the loop field of a 
(collective) field theory of the matrix model (more accurately, is a negative mode part 
of a Bose field) . The field theory is proved to be successful approach to Hermitian matrix 
ensembles [22] • Elsewhere we will develop this approach for the non- Hermitian ensembles. 
The potential is harmonic outside the support of eigenvalues except at infinity where it 
behaves as 2Nhlog \ z\. The current is holomorphic outside the support of eigenvalues. 

In order to obtain the correlations of the potentials, one has to vary the partition 
function by W{z) W{z) + e\og\z — CP where C is a parameter. We denote this 
particular deformation of the potential by 6^: 

6^W{z) =e\og\z-C\\ (16) 

Under this variation the correlation function changes by insertion of the field <f{C)'- 

/z\(A)=e(A^(C))_. (17) 

This is the linear responce relation used in the Coulomb gas theory [TH] . 

While varying the potential it is important to distinguish a harmonic variation of the 
potential W{z, z) —>■ W{z, z) + V{z) + V{z), where V{z) is a holomorphic function. This 
variation does not change AW . To implement a harmonic variation one may extend 
the potential by adding a harmonic function W ^ W + 2TZe Y.k>i ikZ^ and apply the 
operator D{z) = I]fc>i \~^tfc used in Refrs. |Z]-[ini- Then correlators of holomorphic 
parts of the potential 0(2;) = ?iti\og{z — M) are 

(0(2)) = hN log z-n^D{z) log Zn, 

/ n \ n (18) 

( n ) = (-1)"^'" n D{z,) log Z^, n>2. 

\i=l I conn «=1 

Dirichlet boundary problem. We list some elements of the external Dirichlet bound- 
ary problem, which are extensively used below. More details can be found in [221 121] • 

Let Z) be a closed connected domain of the complex plane bounded by a smooth 
curve. Given a real analytic function f{z) in a vicinity of the boundary, we may restrict 
it to the boundary of D. The external Dirichlet problem is to find a harmonic function in 
the exterior of whose value on the boundary is f{z). We call this harmonic function 
the harmonic extension of f{z) to the exterior and denote it by f^{z). It is given by the 

• Dirichlet formula: 

f"iz) = -^fd^,G{z,z')fiz')\dzl 
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where dn' is the normal derivative at the boundary with respect to the second 
variable, with the normal vector being directed to the exterior of the domain D. 

The Green function G{z, z') for the exterior problem is a harmonic function ev- 
erywhere outside D including infinity except the point z = z', where it has a 
logarithmic singularity: G{z,z') log |-2 — z'\ as z ^ z' . If z' oo, then 
(^(z, oo) — log|2;|. The Green function is symmetric in z^z' and vanishes on 
the boundary. In particular, the harmonic extension of log 1-2 — CI in 2; is 

(log \z = log k - CI - G{z, C) + 00). (19) 

If the point C happens to be inside, G(-2, C) in this formula is understood to be null. 

The Green function can be expressed through a conformal map wi^z) from Qj\D 
(the exterior of D) onto the exterior of the unit disk: 

w{z) — w{z') 



G{z,z') = log 



1 — w(z)w(z' 



(20) 



This formula does not depend on the normalization of the map. It is convenient 
to fix w{z) by the condition that it sends infinity to infinity and the coefficient in 
front of the leading term as z ^ 00 is real positive. 

Neumann external jump is the difference between normal derivatives of a smooth 
function at the boundary and its harmonic extension. The Neumann external jump 
operator R acts as follows: 

f{z) ^ {Rf){z) = 9„-(/(^)-/^(^)), z e dD. (21) 

The upperscript indicates that the derivative is taken in the exterior of the bound- 
ary. As it follows from the Dirichlet formula, the Neumann jump is an integral 
operator on the boundary curve with the kernel given by normal derivative of the 
Green function in both arguments: (R/)(z) = dnf{z) + ^ § dndn'G{z, z')f{z')\dz'\. 
In fact this integral is not yet well-defined since the kernel has a second order pole on 
the contour. The operations of taking normal derivative and contour integration do 
not commute. The above formula has to be understood as dn § dn'G{z, z')f{z')\dz'\ 
Alternatively the formula is understood as the principal value integral: 

(R/)(z) = dJiz) + ^VN.jdndn>G{z,z'){f{z) - f{z'))\dz'\. (22) 

Hadamard formula describes deformation of the Dirichlet Green function under 
deformation of the domain. 

A change of the boundary can be characterized by its normal displacement 5n{z), 
such that 5n{z) is a continuous function on the boundary (see Fig. H]). The 
Hadamard formula expresses the deformation of the Green function through 
the Green function itself: 

SGiz,,Z2) = ^f dnG{z,,0dnG{z2,0M0m- (23) 

ZTT JdD 

A corollary of the Hadamard formula is the variation of the boundary value of a 
harmonic function under variation of the boundary. It reads 

6f''iz) = iRf)iz)5niz), z e dD. (24) 
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Figure 1: The domain D and vectors of normal displacement. 

3 Large N limit 

The large N limit is understood as N oo, h while ^A^ is kept finite. The expansion 
in A^~^ is then equivalent to the expansion in h. We call it a semiclassical limit. 

Semiclassical density and the support of the eigenvalues. To elaborate the semi- 
classical limit, we have to find the maximum of the integrand |\l'7v(-Zi, . . . ,-ZAr)P in ©• 
At finite N, it is given by the conditions c?log |\E'iv|/f?^j = for every i: 

N 2 

j=l,^i 

or by the equation 

(v9o(z) + W{z)) = d, (<^o(^) + W{z)) = , zeD (25) 

in the limit. Here ^po{z) = {f{z)) = / log {z—z'l"^ {pi.^)) {z')dPz' and {p{z)) are expectation 
values of the potential and the density in the leading semiclassical approximation. The 
last equation holds inside the support of eigenvalues, D. It does not hold outside it. 
In the semiclassical approximation, D is a domain with a well-defined sharp boundary 
determined by the potential W. An assumption that the support is a compact connected 
domain bounded by a Jordan curve implies restrictions on the potential. The restrictions 
do not reduce the number of parameters in the potential but rather ranges of their 
variation. We do not discuss them here. The assumption is valid, for example, for small 
perturbations of the Gaussian potential. 

The solution for density is obtained by applying to eq. 

r a{z) if zeD 
ipiz)) = (26) 
[ a zeC\D. 

The function a is introduced in ((71). Consequently, (po(z) = / \og\z — z'\'^a(z')(fz' . This 

JD 

function is harmonic in the exterior of the domain except at infinity where it has a 
logarithmic singularity. Eq. (|^ means that inside D, including the boundary, it is equal 
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to —W plus a constant. Since W{0) = 0, the constant is (po{0) = / log and 

JD 

so foiz) + W{z) = (fo{^) for z & D. Moreover, according to (|^. 

d^{^o{z) + W{z)) = Q, zedD, (27) 

so both tangential and normal derivatives of (fo + W at the boundary vanish. 
Since (po{z) — Nhlog |zp is harmonic outside, the harmonic extension of W is 

ly^(z) = ^o(O) - M^) + Nh (log \z\' - (log \z\')"^ , (28) 

whence it follows that 

R(w{z) + Nhlog Izl"^) =0. (29) 

This condition means, in other words, that the domain D is such that the function dzW 
on its boundary is the boundary value of an analytic function in C \ D. The latter, 
together with the normalization condition AWd'^z = —AttHN, determines the shape 
of the support of eigenvalues. 

For the potential of the form (jH} with V{z) = J2k>i tk^^ this condition is somewhat 
more explicit. In this case, the shape of the domain is determined, by the relations 

/ z~''dz = tk, - f d^z = hN. 
Ilk JD vr JD 

such thatvr^A^ is the area of D and —Tiktk are harmonic moments of the domain com- 
plementary to D. The problem is thus equivalent to the inverse problem of 2D potential 
theory. 

Using the relation W{z) = ^po{0) — 'fo{z), it is easy to find the value of /i^ log | \l'iv P 
at the saddle point, which we denote by Fq: 

' \{z)o{z!)d^zd^z' . (30) 



/ / log 


1 


1 






IDJD 


Z 





The leading asymptote of the partition function as ^ — is therefore Zjy — e^'^l'^ . By 
inspection one can check that the average density (j26j) can be also obtained by variation 
of (inni): {p{z)) = 6Fo/6W{z). 

The function Fo plays an important role in the complex analysis. It generates con- 
formal maps from the exterior of D onto the unit disk and gives a formal solution to the 
Dirichlet boundary problem. See Sec. 4 of (TUI for details. 



Variation of the support of eigenvalues. Let us examine the change of the support 
of eigenvalues D under a small change of the potential at fixed A^. In general, we can 
write: J^j^ f{z)d'^z = §qjj 6n{z) f {z)\dz\ for any function /, where 6D stands for a strip 
between the domains D{W + 6W) and D{W) and 6n has the same meaning as in the 
previous section (Fig. [T]). 

The variation of the saddle point condition is most conveniently found from (j27j) : 
S{dn{W + <^o)) = Sn dliW + <^o) + dJ{W + <^o) = 0. 
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where z is on the boundary. Since both tangential and normal derivatives of + (^o 
vanish on the boundary, one writes d\iyV + (/^o) = A(P1/ + y^o) = = —Atkt. This 
gives the variation of the boundary: 

Here RJiy = dn{SW — (SW)^) is the external Neumann jump operator defined by ()2ip . 
(We used the fact that —Sipo is, up to a constant, the harmonic extension of 6W to 
the exterior domain, as is seen from ()28j) .) Note that if the variation of the potential is 
harmonic outside D including infinity the domain does not change. 

Let us check that our result meets the requirement that N stays constant under a 
variation of the potential. Since 



i Sn(z)a(z)\dz\ = — (fdnSW(z)\dz\ = — [ ASW(z)d 
J 4:71 J An Jd 

it is easy to see that 

6 J {p{z))d^z = j 5n{z)(y{z)\dz\+ j^5(y{z)d^z = Q, 

so that is kept constant under this variation. 

Consider special variations of the potential of the form 5cW = e log |z — Then one 
is able to express R(5^iy through the Dirichlet Green function (cf. p9p ). It is convenient 
to introduce the modified Green function: 

giz, C) = Giz, - G{z, oo) - G(cx), C) . (32) 

It is easy to see that 

(The last term in Q vanishes under the normal derivative. It is included for the symmetry 
and future convenience.) 



Variation of the boundary under variation of the size of the matrix. If one 

varies the size of the matrix, N ^ N + 6N, keeping the potential fixed, the support of 
eigenvalues also changes its shape. To find how it grows, we note that 6W = means 
that Sipo{z) = Sipo{0) for all z E D. Plugging here the integral representation of (/Jq, we 
conclude that / log \z~^ — S,^^\5n{^) a{^)\d^\ = for all z in D, with the variation of the 
normalization condition being h6N = § 6n{z)(y{z)\dz\. These conditions are met with 

Sn{z) = - J^^^ dnG{z,oo), (33) 
2na{z) 

where G{z, oo) = — log |u^(-2)|. This describes the interface dynamics in Laplacian growth 
models, known as Darcy's law (cf. |6j). 
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4 Connected correlation functions in the first non- 
vanishing order in h. 



Connected two-point function. In order to obtain the two-point function one has 
to vary the one point function {p{z)) = Q{z; D) o"(z), where 9(z; D) is the characteristic 
function of the domain D: Q{z;D) = 1 for z & D and Q{z; D) = for z ^ D. The 
variation reads 

6{p) = 5aQ{z;D) + a5Q{z;D). (34) 

The second term is locahzed on the contour. Let 6{z; dD) be a ^-function located on the 
boundary of D, defined by the condition that / f{z)6{z;dD)(fz = f{z)\dz\ for any 
smooth function /. It is clear that 6Q{z;D) = 6n{z)6{z; dD). Using the relation 
between a variation of the potential and deformation of the domain, we write: 

And (p) = -A6W e{z; D) + (R 5W) 5{z; dD). 

To keep track of the singular boundary terms, it is helpful to integrate the variation 
of density with some reasonable function / on the plane. This is equivalent to calculating 
5 (tr f{M)) instead of 6 (p). Setting 6W = eg we have: 

— [6{p) fd'z = - [ Ag fd'z + / fiRg) \dz\ . 
e J Jd JdD 

The result is symmetric with respect to f ^ g. With a help of the Green formula, it can 
be expressed through the Bergmann kernel ()22|1 : 

h-'{tTfiM)tTg{M))^^^^ = ^Ue jjj{z)d,g{z)d^z + 

(35) 

+ ^ // f{z)dndn'G{z, z')g{z')\dz\\dz'\. 

Choosing /(z) = log|2;i — 2;^, g{z) = \0g\z2 — zl"^ , we find the pair correlation function 
of the Bose field ip{z): 

^ (<^(Zl)V2(2;2))conn = ^i^U Z2) " log ^^^7^, (36) 

where Q is introduced in (jSH) and logr = lim2^oo(log \z\ + G{z, oo)) is Robin's constant 
(r is the exterior conformal radius of the domain D, see e.g. [21] )• The function in the 
r.h.s. is harmonic outside the domain. If one of the points, say zi, is located inside, one 
sets the corresponding Green functions in (j32j) to be zero. In particular, if both points 
are inside, the correlation function is just 

^ (</^(2l)(^(22))conn = " ^^7^" (37) 

This result is valid for well separated points (|-2i — 2;2p >> NTi/a). 

Taking holomorphic or antiholomorphic derivatives of (jHUj) . we find pair correlations 
of currents: 

( J(^l) J(^2))eo„n = + '^9-.A.G{z,, ^2), (38) 

{Zi — Z2} 
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h'^ {j{z,) J{Z2)) = -7rS^^\z, - Z2) + 2d,,d,,G{zi, z^) (39) 

\ / conn 

(here it is implied that both points are outside). These resuhs resemble the two-point 
functions of the Hermitian 2-matrix model found in [2]; they were also obtained in [T3] 
in the study of thermal fluctuations of a confined 2D Coulomb gas. 

Outside the domain formulas (jHUj) . (jHHj) describe correlations also at merging points 
away from the boundary. In particular, the mean square fluctuation of the current is 

n-Uj\z)) =\{w-z}, zeC\D. (40) 

\ / conn 

Here 

is the Schwarzian derivative of the conformal map w{z). 

These formulas show that there are local correlations in the bulk as well as strong 
long range correlations at the edge. (See for a similar result in the context of classical 
Coulomb systems). At the same time further variation of the pair density correlation 
function suggests that, starting from n = 3, the connected n-point density correlations 
vanish in the bulk in all orders of h (in fact they are exponential in 1/Ti). All the leading 
contribution (of order ^^"'~^) comes from the boundary. 

Note that the result (jH^ is universal in the sense that it depends on the potential only 
through the form of the domain D. The universality holds for any connected correlation 
function of two traces. 

The kernel in the boundary term in (|35|) is the absolute value of the Bergman kernel 
[2H] of the domain C\D aX the boundary. Presumably, this result can be generalized to 
the more complicated case of non-connected supports of eigenvalues, with the boundary 
term being expressed through the Bergman kernel of the Schottky double of the Riemann 
surface C\D. For the Hermitian one-matrix model, where the support of eigenvalues 
shrinks to a number of cuts on the real axis, a similar result was recently obtained in 



Connected three-point function. The 3-point function can be obtained by further 
varying eq. (j35p . Let us first transform the r.h.s., to bring it to the form convenient for 
the varying. Using the Green formula, one rewrites the two point function as an integral 
over the entire complex plane plus the integral over the exterior of the domain 

ATin-^iif{M)iig{M))^^^^ = - f fAgdh+ f {g - g'')Afd'z . 

Jc Jc\D 

The variation of the first term is zero since it does not depend on the contour. The 
variation of the second term consists of two parts: the one coming from variation of 
the boundary and another one from variation of the integrand. The first part actually 
vanishes because the integrand equals zero on the boundary, where g = g^ . The variation 

of the integrand is — / (5g^)Af(fz= f (5g^)(Rf)\dz\ (where we substituted / 

Jc\D J 

f — under the Laplace operator and again used the Green formula). Next, we use 
the Hadamard formula ()23|) to compute the Sg^ , the response of the harmonic extension 
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of the function to a small change of the domain, taken on the boundary of the initial 
domain (j23)- The result is 

47rrt {ti f{M)tTgiM))^^^^ = f{Rf){Rg)5n\dz\. 

Now, let us redefine f = fi, g = f2 and plug (jHH) for the 6n{z) with a function f^: 

5n = Rf'?. Using (IT7I) . we obtain the correlation function of three traces: 

Ana ^ 




The answer is non-universal, i.e., it depends explicitly on the boundary value of the 
Laplacian of the potential. Note that if at least one of the functions fj is harmonic in 
C\D, then the correlation function vanishes (in this leading order in h of course). 

Alternatively, one can apply the Hadamard formula directly to the 2-point correlation 
function of the potentials (p (jHUj) . We have 

((^(Zl)(^(2;2)<^(2;3))conn = —^^3 {v{zi)^{z2)) ^^^^ = 2 — ^3 [g{zi, Z^) + log r] = 

= — ldr.Q{z^,i)dM^2.i)KAi) m = 

€71 J 

The result agrees with the formula for the third order derivatives of the tau-function 
obtained in within a different approach (also based on the Hadamard variational 
formula) . 

For example, in the case W = — l-zp, the support is a disk of the radius R = VhN. 
The conformal map is simply w{z) = z/R and the above formula gives 

(n ,(A,)) ^ -mR^ne ( ' + II « 21. + II » 31 ) . 

conn 

Finite N formulas for the general n-point correlations reviewed in J3] should be able to 
be used to reclaim the same result. 

Variation of contour integrals. One may proceed in the same way to find the con- 
nected n-point correlation functions. Starting from n = 4, however, one encounters 
technical difficulties. 

When passing from 2-point to 3-point functions we transformed integral over the 
boundary into a bulk integral since the latter is easier to vary. Passing to 4-point func- 
tions, one needs to vary the contour integral in (j42|) which does not seem to be naturally 
representable in a bulk form. We have to learn how to vary contour integrals. Here are 
general rules. 

Consider the contour integral of the general form § F{f{z), dnf{z))ds where ds = \dz\ 
is the line element along the boundary curve and F is any fixed function. Calculating the 
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linear response to the deformation of the contour, one should vary all items in the integral 
independently and add the results. There are four elements to be varied: the support of 
the integral §, the dn, the line element ds and the function /. By variation of the / we 
mean integration of the old function over the new contour. This gives § 5ndnF ds. The 
change of the slope of the normal vector results in 5-^ = —ds{5n)^. The rescaling of the 
line element gives 6ds = nSnds, where k{z) is the local curvature of the boundary curve. 
The curvature is k = dO/ds where 9 is the angle between the outward pointing normal 
vector to the curve and the x-axis. The formula n{z) = 9„log \w{z)/w'{z)\ expresses the 
curvature through the conformal map. This formula is useful in calculating the 5-point 
function. We do not attempt to do this here. Another element is the Laplace operator on 
the boundary in terms of normal (dn) and tangential (dg) derivatives: A = d'^ + d] + ndn- 

At last, we have to vary the function / if it explicitly depends on the contour. In 
particular, if this function is the harmonic extension of a contour-independent function 
on the plane, its variation on the boundary is given by (j^^ . 



Connected four-point function. Let us elaborate the case n = 4. In accordance 
with (fTTjl . we find ^n^=i tr /^(M)^ from the response of the r.h.s. of eq. (j^^ to the 
variation of the potential 5W = ef^, so that 



6a = — — A/a , 6n 

471 



Ana 



R/4. 



(43) 



To get the result, we apply the above rules to the contour integral (j42|) . Note that in 
this particular case a change of the slope of the normal vector gives no contribution since 
the functions under the normal derivative do not change along the boundary. Along this 
way we obtain the following result: 

^(ntr/,(M)) =/^EA/. n R/^- 

\j=l ^ i=l k=l,^i 



/\dz\ 
^(9„loga + 2/€)nR// 
^ fe=i 



(44) 



— f \dz\ 



H 



a 



a 



+ [1 ^ 3] + [2 ^ 3] 



The first two terms are explicitly symmetric with respect to all permutations of (1, 2, 3, 4). 
The third one is seemingly not, but in fact it is, as is clear from the Green formula. 

In the 4-point correlation function the first term in the r.h.s. of ()44|) vanishes. For 

cr{z) = l/vr, we get 



n ^(A.) 



\j=i 



\dz\dn f\dz'\ dn'G{z,z') X 



X 



[d„gi\l,z)dngi\2,z)dn'gi\3,z')dn,gi\,,z') + [1 ^ 3] + [2 ^ 3] 

~27r J l^^l'^(^) n dnG{Xj, z), 



(45) 
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where Q is defined by (fH^ and Aj's are assumed to be outside. 



5 Discussion 



We have shown that large scale properties of the normal matrix ensemble are obtained 
from the analytical properties of a curve on the complex plane. The curve bounds a semi- 
classical support of eigenvalues and is determined by the relation dW{z) = —h tij^y 

Large scale correlation functions of the ensemble are objects of the Dirichlet boundary 
problem for the non-compact exterior domain complementary to the compact domain 
D. The two point function is essentially the Dirichlet Green function, the higher order 
functions are related to the deformation of the Green function under deformations of the 
curve. They are determined by successive applications of the Hadamard formula and are 
therefore expressed through the Neumann jump on the curve and through the Bergman 
kernel. 

We expect that other objects of matrix ensembles, such as the genus (1/A^) expansion 
of the partition function are recorded in the analytic properties of the curve. In particular, 
we expect that Fi - a. genus 1 correction to the partition function log (a correction 
of order h!^) is related to the determinant of the Laplace operator in the exterior domain 
C\D. 
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